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In this work we analyze the instabihty conditions for spin-density- waves (SDW) formation on 
the surface of strong topological insulators. We find that for a certain range of Fermi-energies and 
strength of interactions the SDW state is favored compared to the unmagnetized and the uniform- 
magnetization states. We also find that the SDW are of spiral nature and for a certain range of 
parameters a Skyrmion-lattice may form on the surface. We show that this phase may have a non 
trivial Chern- number even in the absence of an external magnetic field. 
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The gapless surface states of strong topological insula- 
tors have drawn a great deal of attention over the past 
few years. These metallic states arise from a topologi- 
cally non-trivial band structure [TJ [2] and are protected 
as long as time reversal symmetry is maintained. This 
topological feature may lead to novel quantum states on 
the surface of these materials, such as Majorana fermions 
and fractional excitations |3]-E]- 

In a previous work [7] it was shown that for a strong 
enough electron-electron interaction, both of the short 
range and the Coulomb types, the surface of a strong 
topological insulator is unstable to the formation of spon- 
taneous uniform magnetization. Other works [SI |9] sug- 
gested that a spin-density-wave state (SDW) is likely to 
occur due to a strong Fermi-surface nesting. In this 
work we explore the conditions under which the SDW 
state is favored over the unmagnetized and the uniform- 
magnetization states. Moreover, we find that for a cer- 
tain range of parameters, a Skyrmion-lattice may form 
on the surface. We also elaborate on the topological 
properties of this phase and claim that a network of one- 
dimensional chiral channels can be established on the sur- 
face. 

We start by analyzing the spin-susceptibility function. 
The non-interacting Hamiltonian of a 3D topological in- 
sulator surface can be approximated by the Fu model 
0: 



non-interacting spin-susceptibility function is given by: 



H voikxCTy 



(1) 



where k± = kx ^ iky^ a are the Pauli matrices, '^o is the 
electron velocity near the Dirac point and A is the warp- 
ing parameter which originated from the cubic Dressel- 
haus spin-orbit coupling of the bulk. Eq. ([T]) introduces 
an energy scale, = {X~^Vq)^^'^ . The Fermi-surface, 
arising from Eq. ([T]), can be classified to three regions 
P^m2] . A circular region for \Ef\ < 0.55£^*, a hexago- 
nal region for 0.55E* < \Ef\ < 0.9£^* and a 'snowflake' 
region for \Ef\ > 0.9E*. Usually, the 'snowflake' region 
coexists with the bulk states and therefore it is less in- 
teresting. Hence, we assume that the Fermi-energy lies 
in the circular or the hexagonal regions. 

We denote the energies and normalized eigenstates of 
Eq. ([T]) as ek,± and /7k, ±. The stationary (a; = 0) and 
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where np is the Fermi-Dirac function at temperature T 
with Fermi-energy ep andFj^,(k, q) = (/7k,s|cr^|/7k+q,s/). 

The eigenvalues of Xo^Iq^ T = 0), for q = (2/ci?, ), as a 
function of the Fermi-energy are presented in Fig. (1(a) ). 
Clearly, as the Fermi energy increases, one of the eigen- 
values becomes much larger than the other two. The 
maximal eigenvalue as a function of q , for a fixed 
Fermi-energy in the hexagonal region, is presented in 
Fig. ( |l(b) ). The susceptibility is maximized at the nest- 
ing vectors of the Fermi-surface in the hexagonal re- 
gion. By analyzing the susceptibility as a function of the 
temperature (see supplementary), we find that the sus- 
ceptibility increases dramatically as the temperature de- 
creases. However, it remains finite for all temperatures. 
This feature arises from the fact that the Fermi-surface is 
not perfectly nested, otherwise we would expect the sus- 
ceptibility to diverge at zero temperature. The Landau 
free-energy up to a second order, with the magnetization. 
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FIG. 1: (a) The zero temperature eigenvalues of Xo^ for a 
fixed q = (2/cf, 0) as a function of the Fermi-energy, where the 
Fermi-energy is measured in terms of ^* = [\~^VqY^'^ and 
the susceptibility is measured in terms of (Avo)~^^^. As the 
Fermi-energy increases, one of the eigenvalues becomes much 
larger than the other two. (b) The maximal eigenvalue of 
as a function of q, for a fixed Fermi-energy {Bp — 0.7 E*), 
where q is measured in terms of (A~^i'o)^^^. The susceptibil- 
ity is maximized at the nesting vectors of the Fermi-surface. 
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niq, as an order parameter, is given by: 

^i = E"^-vr"^q (3) 

q 

Without interactions, as shown in Fig. ( |l(a)| ), the eigen- 
values of the susceptibihty are non- negative. Hence, 
Eq. (|3| is minimized by oiq = for every q. However, in 
the presence of a contact interaction, the susceptibihty 
can be expressed in terms of the non-interacting suscep- 
tibihty using the random phase approximation (RPA), 
X^^ = Xo^[(l ~ 9Xo)~^]^^ ^ where g is the contact- 
interaction strength. Since the non-interacting suscep- 
tibihty is finite for ah temperatures, we conclude that 
a critical interaction at which one of the eigenvalues of 
the susceptibility becomes negative exists for all temper- 
atures. Above that critical value a second order phase 
transition to a SDW ground state occurs. The zero- 
temperature critical interaction for a formation of a SDW 
ground state as a function of the Fermi-energy, as calcu- 
lated in the RPA, appears in Fig. The critical in- 
teraction (mean-field) to have a uniform magnetization 
ground state as calculated in [7] also appears in Fig. ([2|. 
At the circular region, where the Fermi-surface nesting 
is negligible, the uniform state is favored. As the Fermi- 
energy increases to the hexagonal range, the nesting be- 
comes stronger and the SDW state is favored compared 
to the uniform one. We also find that the values of the 
critical interaction for a formation of a SDW ground state 
are quite sensitive to temperature. The zero temperature 
values in Fig. Q are doubled at T ^ £:*/10. 

Above the critical interaction a non-uniform magneti- 
zation will occur on the surface. By denoting the normal- 
ized eigenvectors of the susceptibility as nq, the magne- 
tization becomes 

m(r) (X ^ Re(nq) cos (q • r + (/9q)+Im(nq) sin (q • r + (f^) 
q 

Since all the eigenvectors of the susceptibility contain 
both real and imaginary parts, each q contributes to the 
magnetization a spiral SDW propagating in the q di- 
rection. The magnitude of the SDW is almost constant 
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FIG. 2: The critical interaction to have a SDW/uniform- 
magnetization ground state as a function of the Fermi-energy 
at zero temperature. For a certain range of energies the SDW 
instability is favored compared to the uniform one. 



and its direction varies continuously in space. The rel- 
ative phase between each pair of SDW, Acp^^ is physi- 
cal. Only two of these relative phases are independent, 
and choosing them is equivalent to choosing the origin of 
the SDW. This freedom will be manifested as Goldstone 
modes (phonons). 

Near the phase transition, only the nesting vectors 
{Q^} will contribute to the magnetization. However, 
these SDW do not necessarily coexist. The Landau free- 
energy up to a fourth order is: 

Fl^YI ^Q. I' + I' + ^I'^Q, I Vp,^, I' (4) 

i 

where U and V are four-spin correlation functions. 
Clearly, by tuning the ratio between U and the sys- 
tem can change its phase from a phase where all the spin- 
spirals coexist, to a stripes phase where only one of the 
spin-spiral has a magnitude that is different from zero. In 
general, U and V depend on the dimensionless parameter 
X = g{Xvo)~^^'^ and on the dimensionless Fermi-energy 
(Ef/E^). Though the exact calculation of U and V is 
difficult, a numerical mean-field calculation suggests that 
both phases may occur for different combinations of these 
parameters. 

To the end of this paper, we will focus on the situation 
where the ratio between U and V is such that the favored 
ground state is the coexisting phase. The superposition 
of the non-coplanar spin-spirals form then a triangular 
Skyrmion-lattice on the surface [13[ [14] . 

The components of the magnetization form a periodic 
structure with a triangular symmetry. Each unit cell 
of the triangular lattice holds a Skyrmion with a unit 
Pontryagin-number [15 , Qtop = ±1, and a zero net mag- 
netization. In the core of each Skyrmion the magneti- 
zation points in the direction normal to the surface and 
it rotates continuously towards the opposite direction at 
the edges of the unit cell. An illustration of the mag- 
netization direction in a unit cell appears in Fig. (|3|. 
Clearly, the system chooses spontaneously one of the two 
time-reversal-breaking states which are characterized by 
opposite Pontryagin-number, Qtop = ±1- 

In the u = 1 quantum Hall state [l6l[T? , the Skyrmion- 
lattice phase breaks both the translation and the U{1) 
symmetry of the single-Skyrmion state. Unlike the quan- 
tum Hall case, in this case there is no rotational sym- 
metry and therefore we expect that the only massless 
Goldstone modes will be the standard phonons of a two 
dimensional triangular lattice. 

Motivated by the RPA results, we now discuss the 
mean-field (MF) theory of the Skyrmion-lattice phase. 
This mean-field treatment will allow us a better under- 
standing of the Skyrmion-lattice phase and its topologi- 
cal features. We consider Eq. ([T]) as the non-interacting 
part of the Hamiltonian. We also assume a non- vanishing 
expectation value of the magnetization: 

{^t(r)^^(r))=m(r) = ^m<,e*Q- 

Q 
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where m(r) describes the Skyrmion-lattice that we found 
in the RPA, and Q are the nesting vectors of the Fermi- 
surface in the hexagonal region. In the mean-field ap- 
proximation, the interaction part of the Hamiltonian be- 
comes: 

^x^^ =9Y.\niJ-2gJ2ci{in^-cT) Ck.q (5) 

Q Q,k 

where Ck are the Fourier components of ?/^(r) and g is 
the interaction strength of a contact interaction. The 
full mean-field Hamiltonian, including the kinetic part 
Eq. ([l]), can be solved now by using Bloch theorem, to 
yield a band structure. To the end of this paper, we 
restrict the Fermi-energy to be at the hexagonal range 
and we assume that the magnetization energy is much 
smaller than that scale. The two continuum branches 
of the Dirac-cone split into a periodic band structure 
Fig. Qa)). We define the ground-state MF energy: 

E^f = J2gml+Y^e^,u (6) 

Q n,k 

where the summation of k is over the first Brillouin zone 
and the summation of n is over the occupied bands. By 
minimizing Eq. (|6| with respect to the magnetization, we 
find that a critical interaction exists, above which there 
is a non-trivial minimum to Eq. (|6| . This critical interac- 
tion as a function of the Fermi-energy appears in Fig. ([2|. 
We notice that the MF values of the critical interaction 
agree quite well with the RPA results. We also find the 
phase transition to the gapped phase to be a second or- 
der phase transition where the gap grows as ^yg — As 
shown in [7 , in the known topological insulators such 
as Bi2Te3 and Bi2Se3, these values of critical interaction 
strength are too strong for an instability, albeit not by a 
large factor, making the instability an issue that may be 
relevant for other topological insulators. 

The spectrum of the surface splits into many bands 
which, in general, cross each other. However, we found 
that the two bands near the Dirac-point cross each other 




FIG. 3: The direction of the magnetization at selected points 
in the unit cell of the triangular lattice for the Qtop — 1 phase. 
In the core, the magnetization is in the direction normal to 
the surface. The red line denotes all the points where the 
magnetization is in the plane. 




FIG. 4: (a) An illustration of the transition from a contin- 
uum Dirac cone to a periodic band structure, (b) Magnifi- 
cation of the surface band structure near the Dirac-point for 
Ef — 0.7 E* and g = ^-^gc, where gc is the critical interac- 
tion strength. Where kx and ky are in the first Brillouin zone 
and measured in terms of {X~^voy^'^ and the energy is mea- 
sured in terms of E* = (A^^t'o)^^^. A gap is opened near the 
Fermi-energy at every point in k-space. 



only at the Dirac-point and are gapped from all the 
other bands. The two bands near the Dirac-point for 
representative values of the parameters are presented in 
Fig. Qb)). Clearly, while near the Dirac-point the origi- 
nal branches of the Dirac-cone are almost unchanged, at 
the Fermi-energy the change is dramatic. A gap is opened 
at the Fermi-energy and the surface spectrum becomes 
fully gapped. Moreover, the surface spectrum is not sym- 
metric with respect to the Dirac point, since the mag- 
netization has a non-zero Pontryagin-number (winding) 
which differentiates between states with different chiral- 
ity. 

We now turn to explore the topological properties of 
the surface. In the absence of interactions, the spectrum 
of the surface is gapless. By adding a mass term, which is 
equivalent to a uniform magnetization, a gap is opened 
at the Dirac-point and a half-integer Hall conductivity 
arises [18 . For an infinite Dirac-cone, this half-integer 
Hall conductivity is precisely quantized. However, in real 
systems the Dirac bands have an energy cutoff and there- 
fore the half-integer Hall conductivity is not precisely 
quantized. It deviates from the half-integer value by an 
amount which scales as the gap energy divided by the 
cutoff energy. Moreover, the gap at the Dirac-point can 
be opened due to interactions or by applying an external 
magnetic field. In both cases the gap that is opened is 
rather small (~ ^mev) which makes the observation of 
the half-integer Hall conductivity a hard task, since ex- 
perimentally, it is required to fine tune the Fermi-energy 
to the narrow gap. 

As opposed to this narrow gap, the hexagonal range is 
covering a wide energy window, both in the upper and 
in the lower branch. In the well known topological in- 
sulators this energy window, AE^ex ^ 0.35£^*, is of the 
order of ~ lOO^gy- As shown in Fig. ([2|, by tuning the 
Fermi-energy to any value in the hexagonal range and 
for a strong enough interaction, a gap is opened at the 
Fermi-energy and a quantum Hall effect may be observed. 
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The spectrum of the Skyrmion-lattice phase is always 
gapped at the Fermi-energy and the lattice symmetry of 
this phase yields a Brillouin zone and a wave-function 
of the form of a Bloch wave-function. This allows us to 
define a Chern-number associated with each band: 



A C=±(/7+l) 



(7) 



where, Pj is the projector to the y'^ band. We denote 
the two bands near the Dirac-point as C and where C 
is the lower band and U is the upper band (see Fig. Q). 
Of course, there are many more bands above and below 
these two, both surface and bulk bands. However, the 
C^U bands are gapped from all other bands and we will 
treat them as an isolated band inside the bulk gap. Us- 
me: Eq. 0, we find that when the Fermi-energy is at 
the hexagonal range the Chern-number of this band is 
Cc-\-u = Qtop = ±1 (see supplementary). This Chern- 
number of the combined band is independent of the ex- 
istence of a mass term that introduces an energy gap 
between them. In the presence of a mass term these two 
bands are separated and their total Chern-number is un- 
changed. 

It is instructive to examine how transport on the sur- 
face of a thick slab depends on the Fermi-energy when 
electron-electron interaction is strong enough for sponta- 
neous magnetization to occur. When the Fermi-energy 
is tuned to the hexagonal range in the lower branch, 
a gap is opened beneath the C band and the occu- 
pied bands below the Fermi energy contribute a Chern- 
number ±n. The ± signs refer to the two possible time- 
reversal-breaking states. In general, the contributions to 
this n come from the three dimensional bands within the 
bulk of the slab and the occupied bands on its surfaces. 
However, since the breaking of time-reversal-symmetry 
is limited to the surface, the bulk bands do not con- 
tribute. As the Fermi-energy is increased, the Chern- 
number does not change as long as the system is within 
the Skyrmion phase, since the energy gap does not close. 
Once the Fermi-energy reaches the region where the mag- 
netization becomes uniform, the gap shifts away from the 
Fermi-energy to the Dirac point, and the Hall conductiv- 
ity ceases to be quantized. It becomes quantized again 
when the Fermi-energy is at the Dirac point, in an en- 
ergy gap. Then, the Hall conductivity of the entire slab 
is quantized to an integer, and that of each surface is 
approximately quantized to a half-integer. Finally, when 
the Fermi-energy is increased further, it passes through a 
gapless region in the circular range and gets to the hexag- 
onal range on the positive side of the spectrum. Again, 
the system chooses spontaneously between the two time- 
reversal-breaking states and the occupied bands will con- 
tribute a Chern-number of ±(n+l). The particular value 
of n is not universal, and would depend on details of the 
band structure. We can, however, generally say that the 
difference in Hall conductivity between the Fermi-energy 
lying in the lower and upper hexagonal ranges is an odd 
integer. These two cases are illustrated in Fig. ([5|. 
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FIG. 5: An illustration of the surface band structure in the 
Skyrmion-lattice phase. The valence and conduction bands 
are composed of bulk and surface bands. The bands C and 
U are separated from the other bands by an energy gap. The 
Chern-numbers of the two scenarios differ by an odd integer. 



The analysis above points out to two types of situa- 
tions where the Hall conductivity on a surface may vary 
in space. First, that may happen when time-reversal is 
broken in opposite ways in different regions. And sec- 
ond, that may happen when the Fermi-energy varies in 
space, due to a variation between local density of elec- 
trons. In both cases, the domain walls must carry one- 
dimensional chiral edge channels. Notably, these chiral 
edge modes occur without an external magnetic field. 
Such one-dimensional edge channels on the surface of 
topological insulators may arise as a consequence of a 
magnetic field, as discussed in [20 . Here we point out 
that in the presence of strong enough interactions, they 
may appear without an external field and without the 
need to fine tune the Fermi-energy. An illustration of 
this scenario is presented in Fig. ( 6(a)|6(b")] ). The starting 
point is a topological insulator where the Fermi-energy 
of all the surfaces is tuned to the hexagonal range at the 
lower branch. By tuning the Fermi-energy of selected 
surfaces to the hexagonal range at the upper branch, one- 
dimensional edge channels will appear. This tuning may 
be done by an external voltage. 

To summarize, we considered here the possibility that 
a non-uniform magnetization spontaneously appears on 
the surface of a three dimensional strong topological in- 
sulator due to interactions between surface electrons. We 
assumed that the Fermi-energy is such that the Fermi- 




FIG. 6: An illustration of the one-dimensional edge channels 
on the surface of a topological insulator, (a) The Fermi-energy 
on surface 1 is tuned to the hexagonal range at the upper 
branch while the Fermi-energies of all the other surfaces are 
tuned to the hexagonal range at the lower branch, (b) Same 
scenario as in (a) only now the Fermi-energies on surfaces 1 
and 2 are tuned to the hexagonal range at the upper branch. 
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surface arising from a 2D effective Dirac Hamiltonian 
is strongly nested. We treated interactions within the 
mean- field and the random phase approximations. We 
found that for a strong enough interaction the surface 
may be unstable to the formation of a Skyrmion-lattice 
phase which is a result of a multiple SDW states. We 
found the transition from the unmagnetized to Skyrmion- 
lattice phase to be of second order as a function of in- 
teraction strength. We also found that as long as the 
Fermi-energy is located at the hexagonal range, the de- 
pendence of the critical interaction strength on Fermi- 
energy is rather mild while the temperature dependence 
is rather strong. We showed that even in the absence 
of an external magnetic field a non-trivial Chern-number 
arises from the mean- field Bloch wave function of the sur- 
face in the Skyrmion-lattice phase and that this Chern- 
number can be controlled by tuning the Fermi-energy. 
Using this feature, we suggested that a network of one- 
dimensional chiral channels can be established on the 
surfaces of a strong topological insulator with a strong 
enough electron-electron interaction. 

The authors thank the US-Israel Binational Science 
Foundation and the Minerva foundation for financial sup- 
port. 



Supplementary material 



Temperature Dependence 



0) and non-interacting spin- 



^''(q,T), is given in Eq. (2). 



The stationary {uj 
susceptibility function, Xo 
When the Fermi-energy is at the hexagonal range and q 
is close to one of nesting vectors, one of the eigenvalues of 
Xq^ becomes much larger than the other two. We denote 
the largest eigenvalue of Xo^ as xo- We found that xo 
is maximized at the nesting vectors of the Fermi-surface. 
Now, we would like the examine the temperature depen- 
dence of Xo- In order to do so, we set q = (2A:f, 0). The 
ratio of Xo(^) and Xo{T = 0) as a function of the tem- 
perature appears in Fig. (S-7), where the temperature 



is measured in terms of Clearly, the susceptibility 
increases as the temperature decreases. However, it re- 
mains finite for all temperatures. The critical interaction 
strength scales as Xo hence, its zero-temperature val- 
ues are doubled at T ?^ E'^/IO. The values of xo(q, ^ 
at the other two nesting vectors are identical to Fig. (S- 
[t]) due to symmetry considerations. For other q vectors, 
the values of xo (q, T) are different but the temperature 
dependence is similar. 



Chern-numbers Calculations 

The magnetization of a Skyrmion-lattice is: 

m(r) = M^Re(nQe^-"+^Q) (S- 
Q 




FIG. S-7: The ratio of the largest eigenvalues of x^''(q,7') 
and x^^(Q5 T = 0) as a function of the temperature for a fixed 
q = (2A:f, 0). The temperature is measured in terms of E* = 
(vo/Xy^'^. The susceptibility increases as the temperature 
decreases. 



where the summation is over the nesting vectors {Q}, 
are the eigenvectors of the spin-susceptibility, M is 
the total amplitude of the Skyrmion-lattice and (p^ is 
the phase of each SDW. The mean-field Hamiltonian of 
a strong topological insulator surface with Eq. (S-8) as 
an expectation value of ('^^(r) aijj{r)) is: 



Q 



J2 ClHoC^ -2gJ2 Ciim^ ■ <T)Ck.Q 

k Q,k 

(S-9) 



where Ck are the Fourier components of V^(r), and 
is the non-interacting part which can be approximated 
by the Fu model. Eq. (S-9) can be solved numerically 
via Bloch theorem. That introduces a band structure 
and a set of eigenvectors, \^|Jn,\^)^ at each point in the 
lattice Brillouin zone. The Chern-number calculation can 
be done numerically by defining the projection matrix, 
^n,k = |V^n,k)(V^n,k|, and using: 

Cn = ^^Tr(Pn,k[Sfc,Pn,k,5fc^Pn,k]) 
k 

In order to perform the calculation, we must set the pa- 
rameters of the problem. In particular, we must choose 
the exact form of the magnetization by setting the am- 
plitude, M, and the phases, (/Pq, in Eq. (S-8). We found 
numerically that for any non-zero value of the amplitude 
M, a gap is opened at the Fermi-energy. Hence, any two 
states with two different positive (or two different nega- 
tive) values of M can be connected continuously without 
closing the gap. Consequently, they must share the same 
Chern-number. Moreover, two states that differ only in 
the sign of M are two time-reversal-partners and there- 
fore their Chern-numbers are equal in magnitude and 
opposite in sign. Hence, we conclude that the Chern- 
number of the system can not depend on the magnitude 
of M ^ 0, and that the sign of M affects only the sign 
of the Chern-number. For that reason, we can choose a 
specific value of M without losing the generality. For the 
actual calculation, we chose Ef = 0.7 E"" and g/gc = 1-2 
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where Qc is the appropriate critical interaction. We set 
M to be positive such that Um/Ep = 0.01, where Um is 
the average magnetic energy per unit-ceh: 




where the integration is over a lattice unit-cell and 
is a unit-cell area. 



The phases, (/?q, determine the origin of the Skyrmion- 
lattice. Clearly, the spectrum and the Chern-numbers 
are independent of these phases. Therefore, without loss 
of generality, we set all the phases to be zero. After the 
determination of Ep^ g and the exact form of m(r), the 
numerical calculation of the Chern-number is a straight- 
forward task. 



[1] L. Fu, C.L. Kane and E.J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[2] J.E. Moore and L. Balents, Phys. Rev. B 75, 121306 

(2007) . 

[3] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 

(2008) . 

[4] M. Levin, F. J. Burnell, M. Koch-Janusz and A. Stern, 

Phys. Rev. B 84, 235145(2011). 
[5] B. Seradjeh, J. E. Moore and M. Franz, 

Phys. Rev. Lett. 103, 066402(2009). 
[6] J. Maciejko, X. L. Qi, A. Karch and S. C. Zhang, 

Phys. Rev. Lett. 105, 246809(2010). 
[7] Y. Baum and A. Stern, Phys. Rev. B 85, 

121105(R)(2012). 
[8] Liang Fu, Phys. Rev. Lett. 103, 266801(2009). 
[9] Jian-Hua Jiang and Si Wu, Phys. Rev. B 83, 

205124(2011). 

[10] Y.L. Chen, J.G. Analytis, J.H. Chu, Z.K. Liu, S.K. Mo, 
X.L. Qi, H.J. Zhang, D.H. Lu, X. Dai, Z. Fang, 
S.C. Zhang, LR. Fisher, Z. Hussain, and Z.X. Shen, Sci- 
ence 325, 178 (2009). 

[11] D. Hsieh, Y. Xia, D. Qian, L. Wray J.H. Dil, F. Meier, 
J. Osterwalder, L. Patthey, J.G. Checkelsky, N.P. Ong, 
A.V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y.S. Hor, 
R.J. Cava, and M.Z. Hasan, Nature 460, 1101 (2009). 

[12] D. Hsieh, Y. Xia, L. Wray D. Qian, A. Pal, J.H. Dil, 



J. Osterwalder, F. Meier, G. Bihlmayer, C.L. Kane, 
Y.S. Hor, R.J. Cava, and M.Z. Hasan, Science 323, 919 
(2009). 

[13] T. H. R. Skyrme, Proc. R. Soc. London A 262, 
233(1961). 

[14] S. L. Sondhi, A. Karlshede, S. A. Kivelson and 
E. H. Rezayi, Phys. Rev. B 47, 16419(1993). 

[15] S. M. Girvin, in Topological Aspects of Low Dimen- 
sional Systems, Les Houches Lecture Notes, edited 
by A. Comtet, T. Jolicoeur, S. Ouvry and F. David 
(Springer-Verlag, Berlin/Les Editions de Physique, Les 
Ulis, 2000) ;S. M. Girvin, cond-mat/9907002 (unpub- 
lished) . 

[16] L. Brey H. A. Fertig, R. Cote and A. H. MacDonald, 

Phys. Rev. Lett. 75, 2562(1995). 
[17] R. Cote, A. H. MacDonald, L. Brey, H. A. Fertig, 

S. M. Girvin and H. T. C. Stoof, Phys. Rev. Lett. 78, 

4825(1997). 

[18] X. L. Qi, T. L. Hughes and S. C. Zhang ,Phys. Rev. B 78, 
195424(2008). 

[19] D. J. Thouless, M. Kohmoto, P. Nightingale and M. den 

Nijs, Phys. Rev. Lett. 49, 405(1982). 
[20] M. Sitte, A. Rosch, E. Altman and L. Fritz , 

Phys. Rev. Lett. 108, 126807(2012). 



